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Abstract 

We propose a new constraint(l) corresponding to the FESR (with the moment n=-l) 
free from unphysical regions. Using this constraint(l) together with the constraint(2) 
(with the moment n=l), we search for the simultaneous best fit to the data points 
of a[^t an d ratio up to the SPS energies to determine those values at higher 
energies. We then predict = 107.1 ± 2.6 mb, = 0.127 ± 0.004 at the LHC 
energy(y / i =14 TeV). 
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Recently [1,2], we have searched for the simultaneous best fit of the average of 
pp, pp total cross sections (otot )> an d the ratio of the real to imaginary part 
of the forward scattering amplitude (p^ + ') for 70 GeV < Pi a b < Psps(up to 
the largest momentum of SPS corresponding to y / i=0.9TeV) in terms of high- 
energy parameters constrained by the finite-energy sum rule (FESR) [5] with 
moment n — 1. We then predict and in the LHC (y/s=14TeV) as well 
as high-energy cosmic-ray regions. Block and Halzen[3,4] have also rearched 
the similar conclusions based on duality in a different approach. 



Proposal of a new constraint 

The purpose of this Letter is to propose the other new constraint besides the 
previous one in order to constrain the above parameters. Following ref . [1] , we 
consider the crossing-even forward scattering amplitude denned by 
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F(+ , M = + gM with ImF(+V) = *?po . (1) 



We also assume 



ImF (+) (z/) = ImR(u) + ImF P ,( 



^f C o + Cl log^ + c 2 log 2 ^ + %(-) - (2) 



M 2 \ °M ° MJ M \M 

at high energies for v > N. Here, M is the proton(anti-proton) mass and 
is, k are the incident proton(anti-proton) energy, momentum in the laboratory 
system, respectively. (We use the same notation as Ref. [1,2] in this article.) 
Defining 1 



fM(u) = fM(v)-R(u)-Fp,(u)-FM(0) ~ z/* (0) (a(0) < 0) (3) 
for large value of v, we have obtained [1] in the spirit of P' sum rule [7] 

OP i y 1 7 

ReF«(M) = — y -lmF^ + \u)du + — J a^{k)dk 

o k ^ 

where TV = V^ 2 - M 2 ~ TV. The equation (4) was called FESR(l). 2 This 
FESR suffers from the unphysical regions coming from boson poles below the 
pp threshold. Reliable estimates, however, are difficult. Therefore, we have not 
adopted the FESR(l) in the analysis [1,2]. 

Let us now change our strategy to use FESR(l), not at v — M but at some 
intermediate energy v — v x as 



r, \ 2P A f ulmF^iu) J P N f vk M , . J 

1 M 

1 Although Rei 71 ^" 1 "^ (i/) becomes large for large values of ^, a real constant has to be 
introduced in principle since the dispersion relation for Re.F( + ) (u) requires a single 
subtraction constant F( + )(0)[6,2]. 

2 This FESR(l) corresponds to n = -1[5]. 
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If we choose the value of v\ to be sufficiently large, this constraint is not sensi- 
tive to the unphysical regions (the first term of the right-hand side of Eq. (5)) 
as well as ambiguities from low-energy integrals of pp and pp scatterings. 

Suppose we consider Eq. (5) with N = N x and N = N 2 (N 2 > N x ). Taking 
the difference between these two relations, we obtain 

No No 

- J \ lmRM + m \m) r = ^ J ^/" {v)dv ■ (6) 

Let us call this relation as the constraint (1) which we use in our analysis. This 
constraint gives the relation between high-energy parameters c 2 , Ci, Co, (3p>, and 
the cross-section integrals, and is free from the unphysical regions. 

The general approach 

Besides the constraint (1), we have the FESR corresponding to n — 1[5], 

M N 

J vim F {+ \v)alv + — J k 2 a { t ^{k)alk 
o o 

N N 

= J vim R(v)dv + J vim F P >(v)dv . (7) 
o o 

We call Eq. (7) as the constraint (2) which we also use in our analysis. 

The ReF(+) (is) is calculable from the lm.F^ + \v), Eq. (2) by requiring the 
relation F i+ \-v) = (F^(is))* to hold[8]. Therefore, we obtain [9] 



(+) ReF( + \v 



ImF(+)(i/) 

^(c 1 + 2c 2 log^)-^(^) a5 + FW(0) 



Ka tot 

An 



(8) 



The constraints(l),(2) and the formula of <7 t ^ (Eqs. (1) and (2)) and the 
ratio (Eq. (8)) are our starting points. Among four parameters, c 2 ,Ci,c and 
(3pi, the c and (3p> are represented by the other two (c 2 and ci) parameters 
by two constraints. Then, we search for the simultaneous best fit to a[„t an d 
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pW with three parameters, c 2 ,Ci and F^(0). 

Evaluation of cross-section integrals using experimental data 



The integrals of o-£ t > appearing in RHS(LHS) of Eq. (6)(Eq. (7)) are estimated 
by using experimental total cross sections afot and ofot i n Particle Data Group 
2004[10]. A phenomenological formula, Imf l (u)(= j-crl ot ) = jj(c % + c^log-p + 
41og 2 ^) + g(^)°- 5 + f (^)-o.5 + £(^)-i.5 for i = ^ m is used t0 fit exper _ 
imental <r£otKot) in the region of 2.5 GeV< fc <100 GeV (2.592 GeV< k <100 
GeV). The c 2 and c\ are fixed with the values in our previous analysis, 
(4 , 4) = (0.0479,-0.186) ("analysis 2" in ref.[2]), and the 77(103) points of 
data are fitted with four parameters, c l Q , (3 l , d 1 , f l , respectively. The best-fitted 
values of x 2 (x 2 /(N D -N P )) are 148.5/(77-4) and 71.8/(103-4), respectively for 
pp and pp. The large value of x 2 for pp comes from the inconsistency among 
the data of different experiments. In order to obtain good fit to pp data, we are 
forced to pick up some data points giving large ^-contributions to be removed. 
For this purpose we use statistical method, named Sieve algorithm [11, 3]. In 
this method, by minimizing the Lorentzian squared, Aq = J2^=i l n {l + 7 Ax 2 } 
(not the chi squared, x 2 = EiS A X?) wi th 7 = 0.179[11], a "robust" fit is 
obtained to the same data, where Nd is the number of data points and A% 2 
means the ^-contribution of the i-th point. Points giving Axf > Ax 2 nax in this 
robust fit are regarded as outliers, and removed. We take a cut A^ 2 ^ = 6, 
and seven points are removed. 3 After removing these points, we obtain a 
shifted data set, re-fitted by minimizing the conventional x 2 - As a result 
we obtain renormalized x 2 (i nc l u ding the factor R = 1.140[11]), Xpp/(^D- 
Np)=36. 6/(70 — 4). Finally we obtain the successful fits for both afot an d 
aZ- The (4, fi\d\ f) = (6.34, 11.08, 5.28, -0.15)( (6.32,4.25,-12.64,24.4) ) 
are obtained for i = pp{pp). 

We take the values of parameters appearing in Eq. (6) as (Ni,N 2 ,ki) = 

(10, 70, 40) GeV, where v 1 = \Jk\ + M 2 . The above phenomenological fits give 

the cross-section integrals ^ -J^a\ ot {k)dk = 242.57 ± 1.00 (219.04 ± 

0.47) GeV -1 with for i = pp(pp), where the errors correspond to the one- 
standard deviations. By averaging these values we obtain 



2^(^)^ = 230.81 ±0.55 GeV" 1 . (9) 
i 

iVi 

We can also evaluate the cross-section integral in Eq. (7). We devide the 



3 (A;(GeV),fjf o p t (mb))=(2.5,74.9±1.0), (3.54,69.7+0.5), (3.6,76.2±1.8), (4.,71.± 1.), 
(4.015,66.84± 0.32), (4.3,60.6±0.8), (9.14,57.51+0.73) are removed. 
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region of integral into two parts, ± J* k 2 a\ ot (k)dk = ± J k 2 al ot (k)dk + 
fJ^i k 2 a\ ot (k)dk, and the integral in higher energy-region (the second term) is 
evaluated by using the phenomenological fit in the same manner. The integral 
in lower energy region(the first term) is evaluated by using experimental data 
directly: Each datum is connected with the next point by a straight line in 
order, and the resulting polygonal line graph gives the relevant integral. (The 
details of this procedure are explained in our previous works [1,2].) By taking 
the N as 10 GeV and N P P = 4.7GeV, we obtain ^ ' k 2 aZ(k)dk = (522.22± 
1.91) + (3499.44 ± 14.22) = 4021.66 ± 14.35GeV for i = pp. By taking 

N P P = 4.966GeV, ^ jf ' k 2 a'Z(k)dk = (357.24 ± 0.89) + (2411.26 ± 3.50) = 
2768.50 ± 3.61GeV for i = pp. By averaging them we obtain 



This value is consistent with our previous estimate, 3403 ± 20GeV[l], which is 
evaluated by using the area of the polygonal line graphs up to k—N (=10GeV) . 
In our present estimate, both of the integrals are estimated with small errors 
less than 0.3% . 

FESR as two constraints 

By using the integrals, Eqs. (9) and (10), we obtain the constraints(l) and (2) 

as 

constraint(l) 3.316/3^ + 31.98c + 141. l Cl + 610.9c 2 = 230.81 , (11) 
(normalized 0.104/3 P , + c + 4.41ci + 19.1c 2 = 7.22 ), 
constraint^) U0.7[3 P , + 383.6c + 781.6ci + 1635.c 2 = 3395.1 , (12) 
(normalized 0.367/3 P / + c + 2.04 Cl + 4.26c 2 = 8.85 ), 

where we neglect the errors of cross-section integrals, and regard Eqs. (11) and 
(12) as exact constraints. The equations are also rewritten in the form with 
the coefficient of Co normalized to unity in the parenthesis. Solving these two 
equations, we obtain the constraints for Co and f3pi as 

Co = Co ( C2) Cl ) = 6.574 - 5.348c! - 24.95c 2 , 
p pi = (3 P , ( C2 ,ci) =6.206 + 9.025ci + 56.40c 2 . (13) 

Analysis and result 

In the actual analysis we fit the data of ReF^ + \u) instead of p( + \ We made 




(10) 



o 
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(a) a™ (b) P (+) 

Fig. 1. Predictions for and pW; (a) Total cross section 
versus logio-FLft/GeV, (b) gives the p^ + \— Re F^/Im F^) versus E cm in 
terms of TeV. The fit is done for the data up to the SPS energy, in the region 
70(10)GeV< k < 4.3xl0 5 GeV (11.5(4.54)GeV < ^7s < 0.9TeV) for a t ( +V +) ) 
which is shown by horizontal arrow in each figure. Vertical arrow represents 
the LHC energy ^/s=UTeV, corresponding to fc=1.04x 10 8 GeV. The thin dot- 
dashed lines correspond to the one standard deviation of c 2 , given with the 
parameters (c 2 , c u c , (3 P ,, F(+>(0)) = (0.0464 ± 0.0038, -0.158 q= 0.057, 6.26 ± 
0.21, 7.40 =F 0.30, 10.18 =F 0.27) . 

(Jtot and ReF^ + ) data points by averaging the original data given in ref.[10]. 
The detailed explanations for the treatment of data are given in ref.[l]. There 
are 17 data points of above 70 GeV and 10 data points of ReF^ above 
10 GeV up to SPS energy y^s = 0.9 TeV. They are fitted simultaneously with 
parameters c 2 ,ci for a^t an d c 2 , ci, F( + )(0) for ReF^ + \ 

The results are shown in Fig. 1. The % 2 /d.o.f is 10.80/(27 — 3),which is less 
than unity. The respective x 2 -values devided by the number of data points for 
cr£ ) and are xl/ N * = 5.64/17 and x 2 P / N P = 5.16/10, respectively. The fit 
is successful. 



The values of parameters are given in Table 1. The result is compared with 
the previous analysis ("analysis 1" in ref.[2]), where only Eq. (12) is used as a 
constraint. The values of c 2 have to be noted since the high-energy behaviours 
of <j( + ) and p( + ) are most sensitive to c 2 . The "analysis 1" gives c 2 = 0.0466 ± 
0.0047. The present value of c 2 in Table 1 is consistent with the previous one, 
although the error does not reduce so largely. 

In order to check our result, we take another value of k±, A;i=80GeV for con- 
straint^) with the other parameters to remain unchanged ((iVi, N 2: iV)=(10, 
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Table 1 

Values of parameters in the best fit using the constraint(l) with (Ni, N 2 ,ki) 
= (10,70,40)GeV and the constraint(2) with iV=10GeV. The c and (3 P > are rep- 
resented by C2 and c\ through the constraints, and the fit is performed by using 
three parameters C2, c\ and F^ + \0), of which errors are given by the x 2 function 
X 2 ( c 2) ci, F^ + )(0)). Conversely by solving the constraints for C2 and ci, the x 2 func- 
tion is represented by cq, (3p> and F( + )(0). The errors of cq and (3p> are, thus, 
obtained. 



C2 


ci 


CO 


Pp> 


f(+)(o) 


0.0464 ± 0.0038 


-0.158 ±0.057 


6.26 ±0.21 


7.40 ± 0.31 


10.18±1.70 



70, 10)GeV). In this case the constraint(l) becomes — 3.628/3p/ — 27.68co — 
113.0ci-463.3c2 = -203.93 ±0.34GeV~ 1 (normalized:0.131 ( 5 P /±c ±4.08ci± 
16.7c2 = 7.37±0.01). By using this together with the constraint (2) (Eq. (12)), 
the analysis is done in the same way. The value of c 2 is obtained as c 2 = 
0.0472 ±0.0036 with x 2 /d.o.f = 10.84/(27-3). This c 2 is almost the same as 
in Table 1, and the result is considered to be almost independent of the value 
of hi. 

We have also checked the dependence of the errors of cross-section integrals, 
Eqs. (9) and (10). In the case when a larger value of the integral, 230.81±0.55 
(3395.1±7.4), is used for constraint(l) (constraint (2)) with the other integral 
to remain the same, we obtain that the best-fit value of C2 is 0.0485(0.0461) 
with x 2 = 11.37(10.88). The deviations of c 2 from the original value 0.0464 are 
0.0021 (-0.0003). They are small, compared with the statistical error 0.0038: 
about (21 2 ±(-3) 2 )/38 2 = 30%. So, we can regard Eqs. (11) and (12) as exact 
constraints. 

Special attention has to be paid when the two constraints are employed to 
constrain the values of c and (3 pi. As shown in Eqs. (11) and (12), the 
constraints(l) and (2) take the co-normalized forms 0.104/3p/ + cq + • • • and 
0.367/3p/ ± Co ± • • •, respectively. They are linearly independent and we have 
obtained the meaningful result. If the parameters are badly taken so that two 
constraints are not sufficiently linearly independent, the result becomes mean- 
ingless. For example, in case of A;i=60GeV, the constraint(l) has a normalized 
form 0.314/3p/±co±- • •. It is quite close to the constraint (2), and this selection 
of parameters are not suitable for the analysis. In the present analysis, in order 
to obtain sufficiently independent constraint, we have taken much larger value 
of ]V2(=70GeV) for constraint(l) than iV(=10GeV) for constraint (2). 

It is pointed out [4] that there are strong resemblances between our approach 
and the one by Block and Halzen[3]. They estimated the values of experimental 
even-cross section <7 even (i/ ) and of its derivative ^(JJm) I v o a ^ a certain energy 
v — i/q — 7.59GeV(A/s = 4GeV) by using a local fit. These two quantities 
are used as constraints for parameters C2, ci, c and (3pi, and c and (3 pi are 
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represented by C2 and c±, similarly to our Eq. (13). They have shown in Ref.[4] 
that the constraint for o" cvcn (z/ )) which gives 8.67 = c + 2.091ci + 4.371c 2 + 
0.3516/3p/, 4 is very close to our FESR(2), Eq.(12). Numerical difference seems 
very small at a first look, but this difference is physically very important, since, 
in the former, the constraint is obtained at one point v = u in asymptotic 
region of a[^t , while, in the latter, all the information in low-energy resonance 
region is included in the integral of a[„t taken from k = to 10 GeV. 

By using the values of parameters in Table 1, we can predict the a[„t and 
at Tevatron-collider energy(A/i=1.8TeV) and LHC energy ( v /s=14TeV). 



(tS = 75.82 ± 1.02mb (y/s = 1.8TeV), 107.1 ± 2.6mb (Js = 14TeV) 

(14) 

= 0.136 ± 0.004 (v/i = 1.8TeV), 0.127 ± 0.004 (^ = 14TeV) 

where the relevant energies are very high, and the a[^t an d p^ can be regarded 
to be equal to the a PP t and p pp . 

Our predicted values are almost the same as the previous ones[2]. They are 
consistent with the recent prediction by Block and Halzen[3] a PP t — 75.19±0.55 
mb, pP p = 0.139 ± 0.001 at Tevatron energy ^ = 1.8TeV, and af p t = 107.3 ± 
1.2 mb, p pp = 0.132 ± 0.001 at LHC energy Ifs = 14TeV. They also analyzed 
the crossing-odd amplitude and obtained smaller errors compared with ours. 
Our prediction has also to be compared with Cudell et al.[12] a pp t = 111.5 ± 

l-2 sys t -2:l stat mb, p pp = 0.1361 ±0.0015 syst iaoo25 stat , whose fitting techniques 
favour the CDF point at y/s = 1.8TeV. 

Finally we emphasize that our present analysis with two constraints is inde- 
pendent of the previous one [1,2] with one constraint. Although the high-energy 
parameters are strongly constrained by two FESR, Eqs. (11) and (12), in the 
present analysis, the result is almost the same with the previous one [2]. 

It is worthwhile to point out the followings: 

1. Both of the parameters Cq, (3 pi are constrained as Cq = Co(c2,Ci) and 
(3 pi = (3p'(c2,Ci) through FESR (namely duality). 

2. And, then the high-energy behaviours of have been x 2 fitted in terms 
of the parameters C2, c\ since is most sensitive to C2- 
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4 In ref.[4], the constraint is given in unit of mb and the LHS is given as 48.58 mb, 
which is replaced by 8.67 here in our notation where q's are dimensionless. 
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